Nucleosynthesis in Supernovae 



C. J. Horowitz* and Gang Li^' 
Nuclear Theory Center 
2401 Milo B. Sampson Lane 
Bloomington, Indiana 47405 
(February 1, 2008) 

Core collapse supernovae are dominated by energy transport from neutrinos. Therefore, some 
supernova properties could depend on symetries and features of the standard model weak interac- 
tions. The cross section for neutrino capture is larger than that for antineutrino capture by one term 
of order the neutrino energy over the nucleon mass. This reduces the ratio of neutrons to protons in 
the //-driven wind above a protoneutron star by approximately 20 % and may significantly hinder 
' r-process nucleosynthesis. 
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' Core collapse supernovae are perhaps the only present day large systems dominated by the weak interaction. They 
. are so dense that photons and charged particles diffuse very slowly. Therefore energy transport is by neutrinos (and 
' convection). 

. We beleive it may be useful to try and relate some supernova properties to the symmetries and features of the 
On ' standard model weak interaction. Parity violation in a strong magnetic field could lead to an asymmetry of the 
, explosion Indeed, supernovae explode with a dipole asymmetry of order one percent in order to produce the very 
high 'recoil' velocities observed for neutron stars However, calculating the expected asymmetry from P violation 
has proved complicated. Although explicit calculations have yielded somewhat small asymmetries |3|-^ it is still 
possible that more efficient mechanisms will be found, 
t — In this letter we calculate some effects from the difference between neutrino and antineutrino interactions. In 

Quantum Electrodynamics the cross section for e^p is equal to that for e'^p scattering (to lowest order in a). In 
contrast, the standard model has P-nucleon cross sections systematically smaller than ^/-nucleon cross sections. 
However at the low v energies in supernovae, time reversal symmetry limits the difference between v and 9 cross 
, sections. Time reversal can relate v — N elastic scattering and D — N where the nucleon scatters from final momentum 
Pf to initial momentum pi. If the nucleon does not recoil then the v and D cross sections are equal. Thus the difference 
between i/ and i> cross sections are expected to be of recoil order E/AI where E is the neutrino energy and AI the 
nucleon mass. This ratio is relatively small in supernovae. However the coefficient multiplying E/M involves the 
large weak magnetic moment of the nucleon (see below). 

The standard model has larger i' cross sections than those for D. For neutral currents, this leads to a longer mean 
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, free path for ly^ compared to I'x (with x=/i or r). Thus even though i/^ and i^x are produced in pairs, the antineutrinos 
' escape faster leaving the star neutrino rich. The muon and tau number for the protoneutron star in a supernova could 
be of order lO^** Supernovae may be the only known systems with large p and or t number. For charged currents, 
the interaction difference can change the equilibrium ratio of neutrons to protons and may have important implications 
for nucleosynthesis. We discuss this below. To our knowledge, all previous work on nucleosynthesis in supernovae 
assumed equal v and P interactions (aside from the n-p mass difference). 

The neutrino driven wind outside of a protoneutron star is an attractive site for r-process nucleosynthesis Q . Here 
nuclei rapidly capture neutrons from a low density medium to produce heavy elements This requires, as a bare 
minimum, that the initial material have more neutrons than protons. The ratio of neutrons to protons n/p in the 
wind depends on the rates for the two reactions: 

i^e + n ~> p + e~ , (la) 
i>e + P n + . (16) 
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^We expect the difference for charged current interactions to be of the same order if one can neglect the neutron-proton mass 
difference. 
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The standard model cross sections for Eqs. (la, lb) to order E/A4 are, 

with G the Fermi constant (and 9c the Cabbibo angle), E^. — E zL A the energy of the charged lepton and A = 1.293 
MeV is the neutron-proton mass difference. The plus sign is for Eq. (la) and the minus sign for Eq. (lb). We use 
ga ~ 1.26.^ Equation (2) neglects small corrections involving the electron mass and coulomb effects (see for example 
1^) while the finite nuclcon size only enters at order (E/M)"^. 
We refer to the 7 term as a recoil correction. It is the same for i' and D. 

7= (2 + 10g2)/(l + 35^) «3.10 (3) 

Finally, the S term involves the interference of vector (I-I-2F2) and axial (ga) currents. This violates P, which by CP 
invariance also violates C. This increases the v and decreases the 9 cross section. 

4.9,(1 + 2F2)/(l + 3g^)« 4.12 (4) 

Here F2 is the isovector anomalous moment of the nucleon. (This is the weak magnetism contribution.) 
We average Eq. (2) over the Ve spectrum to get, 



<a>, = ^ 1 + 3<?^ <^>e 1 + 2- + ao— 1 + <5 - 7 "277 , 5 

for Eq. (la). Here the mean energy e is defined as, 

t^<E'^>l<E>, (6) 

and a2 is a shape factor a2 = <E^><E> / <E'^>'^ . Finally ao = <E'^> / <E>'^ and <E^> are the ith energy moments 
of the Ve spectrum. Note, e « 1.2<E>. 

Likewise, averaging over the De spectrum for Eq. (lb) gives, 

G^cos'^Oc, 9s - r A A2„ , e, 

<^>P = + 3gl)<E>e[l - 2- + ao^][l -{5 + 7)02^], (7) 

with the mean antineutrino energy e = <E^>/<E> and <E^> the ith moment of the Ue spectrum. We assume 
similar shape factors 02 and ao for 9^ and Ve- The shape factor 02 — 1.23 (1-15) for a Fermi Dirac distribution 
with chemical potential fi = rjT^ and temperature Ty for rj = (3.5). See for example pO| . For simplicity we adopt 
a2 =ao = 1.2 [[ 

The equilibrium electron fraction per baryon Yg (which is equal to the proton fraction assuming charge neutrality) 
is simply related to the rate A for Eq. (lb) divided by the rate A for Eq. (la). 

re = (l + ^)-^ (8a) 



This assumes the neutrino capture rates dominate those for other reactions. Reference |11 contains some discussion 
of the small corrections from capture. The ratio n/p is. 



n 1 

Taking the ratio of Eq. (7) to Eq. (5) gives. 



(86) 



Yc={l + ^QC] . (9) 



^ Note in principle, there is another correction to Eq. (2) from the thermal motion of the nucleons. This is of order T /M and 
increases both the v and u cross sections. However we assume the temperature in the wind T is much less than the neutrino 
sphere temperature and neglect this term. 

^The coefficient ao only makes a very small contribution and our results are insensitive to its value. 
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Here L^^ (Li?^) is the Ve (i^e) luminosity, Q is the correction from the reaction Q value, 



and the C violating term, Eq. (4), contributes the factor C, 

l + (<5-7)«2ir' ^ ^ 

Note, the recoil term 7 makes a small but nonzero contribution to Eq. (11) because the ly and D energies are different. 
Simply evaluating Eq. (11) for typical parameters yields C ~ 0.8. Thus, the difference between v and v interactions 
reduces the equilibrium n/p ratio by approximately 20 %. This is a major result of the present paper and will be 
discussed below. 

Figure 1 shows the values of e and e necessary for = 0.5. We assume equal luminosities Li,^ — Lp^. The region 
to the upper left is neutron rich and to the lower right proton rich. The conditions for ~ 0.5, assuming C = 1 in 
Eq. (9), are indicated by the dotted line. Including C shifts the conditions for Y^ = 0.5 to the solid line. Thus the 
difference in v and 9 interactions converts the region between the solid and dotted lines from neutron rich to proton 
rich. 

We also show in Fig. 1 the values of e and e from a supernova simulation by J.R. Wilson as reported in ref. p2[ . 
The symbols show how the mean energies evolve with time. As the protoneutron star becomes more neutron rich, 
the opacity for De decreases because there are fewer protons. This allows the i>e to escape from deeper inside the hot 
protoneutron star. Therefore e increases with time. Without C the wind starts out with Ye ~ 0.5 and then becomes 
neutron rich. With C the wind starts out proton rich and ends up with Yg w 0.5. If Lp^ w the wind is never 
significantly neutron rich. If Lp^ « I.IL^^ the wind will end slightly neutron rich. However, n/p is still 20 % lower 
with C than without. For example, if Y^ drops as low as 0.42 in a model without C it will only drop to approximately 
0.48 when the difference between v and 9 interactions is included. With this increase in Yg, it is very unlikely that 
successful r-process nucleosynthesis can take place in the wind of this or similar models. 

Note, we are being slightly inconsistent to include the E/M term in Eqs. (2,4) for the neutrino absorption while it 
is not included in the simulation used for e and e. Indeed this term could change the location of the neutrino spheres 
and slightly increase e and decrease e. This could cancel a small part of the effect on the n/p ratio. However, our 
preliminary estimates suggest this change in the spectrum is very small. Including the term in a full simulation would 
be useful jl^] . For completeness we give a C violating term for neutrino-electron scattering NFS which may be useful 
for calculating differences between the Vx and i?x spectrum. 

The total cross section ag for NFS (see ref. for example) is expanded in powers of E/Ep where Ep is the 
electron Fermi energy. To order [E/EpY , 

G'^ E'^ 

with 6f. = AcyCa/3{cf, + c^) and the plus sign is for v and the minus sign for D. The couplings are Cy — 2sin^0i4' ±1/2 
and Ca = ±1/2. Here the plus sign is for Vg and the minus sign for v^- The C violating coefficient w 0.55 for and 
« 0.1 for Although this term is nominally of larger order, E/Ep for NES than E/M for nucleon scattering, the 
coefficient is smaller Se -C S. Therefore we do not expect large differences from NES (except perhaps at low densities). 

With the approximately 20 % reduction in n/p from the difference between 1/ and D interactions, there appears to 
be very serious problems with r-process nucleosynthesis in the wind of present supernova models. In addition to the 
initial lack of neutrons, one has to overcome the effects of neutrino interactions during the assembly of a particles 
and during the r-process itself | [T5[ |. These further limit the available neutrons per seed nucleus. Thus, it is unlikely 
that present wind models will produce a successful r-process. Of course, the wind in supernovae may not be the 
r-process site, although this may be unappealing (see for example |^,|6|). If the wind is not the site, one must look 
for alternative environments. 

However, the effects of neutrino interactions may be very general. The only requirement is that energy transport 
from neutrinos plays some role in helping material out of a deep gravitational well. Given this, it is quite likely that 
the n/p ratio will be determined by the relative rates of Eqs. (la, lb). Therefore differences in v and P interactions 
may be important for just about any nucleosynthesis site that involves neutrinos. Indeed, Haxton et al. p^ ] claim the 
abundance of isotopes produced by neutrino spallation imply significant neutrino fluences during the r-process. 
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If the u-driven wind is the r-process site, it is very hkely, present models of the neutrino radiation in supernovae are 
incomplete. The high values of make it almost impossible to have a successful r-process by only changing matter 
properties, such as the entropy. The neutrino fluxes will (almost assuredly) need to be changed. 

Changes in the astrophysics used in the simulations or new neutrino physics such as neutrino oscillations jist could 
change e, e and or the luminosities and lead to a more neutron rich wind. The oscillations of more energetic Vx 
with Ve could increase e. However, we have some information on the Ve spectrum from SN1987a |^9|. Thus one can 
not increase e without limit. Indeed if anything, the Kamiokande data suggest a lower e. Any model which tries 
to solve r-process nucleosynthesis problems by increasing e should first check consistency with SN1987a observations 
. Alternative modifications could include oscillations of Ve to a sterile neutrino or a lowering of e. (However, we 
know of no model which lowers e.) Whatever the modification of the neutrino fluxes, one will still need to include the 
differences between v and j> interactions in order to accurately calculate n/p. 

In conclusion, supernovae are one of the few large systems dominated by energy transport from weakly interacting 
neutrinos. Therefore, some supernova properties may depend on symmetries and features of the standard model weak 
interactions. The cross secton for neutrino capture is larger than that for antineutrino capture by a term of order 
the neutrino energy over the nucleon mass. This difference between neutrino and antineutrino interactions reduces 
the ratio of neutrons to protons in the i/-drivcn wind above a protoneutron star by approximately 20 % and may 
significantly hinder r-process nucleosynthesis. 

This work was supported in part by DOE grant: DE-FG02-87ER40365. 
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FIG. 1. Mean antincutrino energy e vs mean neutrino energy e, see Eq. (6). The solid line indicates an equilibrium 
electron fraction Y^, — 0.5 including the difference between v and v interactions, C term in Eqs. (9,11), while the 
dotted line shows Yg — 0.5 without this term. The symbols are the mean energies of a simulation by J.R. Wilson as 
reported in ref. for the indicated times in seconds after collapse. The i/-driven wind is neutron rich in the upper 
left of the figure and proton rich in the lower right. The region between the dotted and solid lines is converted from 
neutron rich to proton rich by the C term. 



5 



